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cessing such as quantum teleportation [TJ, superdense coding [2], and quantum key distribution 



The dynamical properties of quantum entanglement in two effective two-level trapped ions 
interacting with a laser field are studied in terms of the negative eigenvalues of the partial 
transposition of the density operator. In contrast to the usual belief that destroying the en- 
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tanglement can be observed due to the environment, it is found that the Stark shift can also 

G ■ 

produce sudden death of entanglement and longdived entanglement between the qubits that 
' are prepared initially in separable states or mixed states. 
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. As a promising resource, quantum entanglement plays a key role in quantum information pro- 
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Qlh, [3J. However, in the real world, quantum information processing will be inevitably affected by 

the decoherence that destroys quantum superposition and quantum entanglement. The extent 
to which decoherence affects quantum entanglement is an interesting problem [I], and many 
researchers study it extensively based on various models [5j El [7J [8] with the point of view of 

> ■ 

■ environment induced decoherence. 

The decay of entanglement cannot be restored by local operations and classical commu- 
nications, that is one of the main obstacles to achieve a quantum computer [9j. Therefore it 
becomes an important subject to study the loss of entanglement [U [TTJ, [T3j [TJJ [15], H5] . Quite 
recently, by using vacuum noise two-qubit, entanglement terminated abruptly in finite time has 
been performed [3] and the entanglement dynamics of a two twodevel atoms model have been 
discussed [TU [T2] . They called the non-smooth finite-time decay entanglement sudden death. 
Although entanglement can be realized in different ways in experiments, how we can preserve 
it is still a big challenge for current technology [17] . Because for open system, entanglement is 
fragile and decays exponentially, it is often thought as similar as quantum decoherence. Most 
of the authors who have treated this problem have dealt with the case in which the Stark shift 
has been ignored [H [10] . However, in reality it cannot be ignored. The main aim of the present 
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paper is try to answer the following question: what happens to two qubits entanglement when 
we consider the Stark shift and different initial state setting in the absence or presence of the 
decoherence? 

We present an explicit connection between the initial state setting, Stark shift and the 
dynamics of the entanglement. We give a condition for the existence of either entanglement 
sudden death or long-lived entanglement. In particular, a quantitative characterization of a 
general system of two three-level trapped ions interacting with a laser field is presented. We 
present various numerical examples in order to monitor the partial transpose of the density 
operator and entanglement dynamics. In principle, by proper adjustment of the initial state 
parameters, we can always find suitable values of Stark shift which can be use to suppress the 
decay of entanglement. 

2 Model 

The physical system on which we focus is two effective two-level harmonically trapped ions 
with their center-of-mass motion quantized. The electronic levels \a) , \b) and |c) are assumed 
to be metastable and coupled via a laser field of the form [T5] 

E(x,t) = E ex.p[i(kx - ojt + (p(t))], (1) 

where Eq is the strength of the electric field, k is the wave vector of the driving laser field, x is 
the position operator associated with the center-of-mass motion and u is the laser frequency. 
We denote by 4>(t) the fluctuations in the laser phase. Therefore, we can express the center-of- 
mass position in terms of the creation and annihilation operators of the one- dimensional trap 
namely x = yhj (2mu s )(a) + a) — Ax(a) + a). We denote by a and a) the annihilation and 
creation operators and u s is the vibrational frequency related to the center-of-mass harmonic 
motion along the direction x. In the absence of the rotating wave approximation, the trapped 
ions Hamiltonian may be written as 

H = H + H 1} (2) 
where H = t\ijj s a)a + hw a \a) (a\ + hwb\b) (b\ + fuv c \c) (c|, and 

2 

H t ^hJ2 (\f e - i{ - k ^- M ^S^ + xfe-^-^^S^ + c.c.) . (3) 
i=i 

We denote by = (b\p^\c) E i (X^ = ( a \p^\c) Eq 2 ) the Rabi frequency characterizing the 
coupling strength (products of dipole matrix elements and amplitudes of the incoming fields), 
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where p {i) is the dipole moment operator. As usual, to describe this system we use the operators 
= \l^)(m^\, (l,m = a, 6, c and i = 1,2). For the sake of simplicity (but without loss of 
generality), we have assumed to deal with the case in which <p l = 0, 2 = an d the level |c) is 
assumed to be dipole-coupled to both the levels \a) and \b) via a far detuned laser field. While 
this is straightforward, it is often the case that it is simpler to work in the interaction picture in 
which the Hamiltonian (j2J) evolves in time according to the interaction with the vacuum field. 
If we express the center-of-mass position in terms of the creation and annihilation operators, 
the interaction part of equation (T5]) becomes 

2 

H int = h*^(sg + S$>) 
i=i 

+h (x^e-W^S® + A^e-^+^W + c.c.) 

+h (xfe-'^^S^e-^ + X^e'^^S^e-^ + c.c.) , (4) 



where r\ = fcy iMj ' * s ^ e Lamb-Dicke parameter and A is the detuning. In equation (py) 
the time-dependent factor is eliminated in the interaction picture, since uo c — (ujb + A) = 
uj s and uj c — (u a + A) = u s (degenerate levels). Making use of the special form of Baker- 
Hausdorff theorem [19] the operator exp[i?7(a^ + a)] may be written as a product of operators 
i.e. exp(i?7(a^ + a)) = exp ^^-[a^,a]) exp [irja^ exp (irja) and assume the Lamb-Dicke regime 
with small rj. In order to obtain this we detune the laser frequency u to the first vibrational red 
sideband. Also, we apply the rotating wave approximation discarding the rapidly oscillating 
terms and selecting the terms that oscillate with minimum frequency [2U]. In these limits we 
can expand the interaction Hamiltonian to lowest order in rj. The resulting Hamiltonian may 
be written as 

2 

H int = hAj2(s2 + SV) + h((?hlsU+(UtfSW+c.c.) 

+H(C?a^t ) e- i UC? ) a'S^ + c.c.), (5) 



with a new coupling parameter Q"^ that includes the Dicke parameter in its definition. The 
analysis of such a Hamiltonian model can be carried out, providing eliminating of the non- 
resonantly coupled atomic level |c) adiabatically, due to the large detuning, the transitions for 
instance from the level \a) to the level |c) are very fast and immediately followed by decays on 
the atomic level Therefore, considering only coarse grained observables, meaning that the 
system is observed at a rough enough time scale, effectively eliminates the far detuned level, 
namely, at such a time scale, the only observables and hence meaningful dynamical behaviors, 
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involve levels \a) and |6) as a result of time averaging second order processes having |c) as 
an intermediate virtual level. This procedure then suppresses the fine dynamics, that is it 
sacrifices any information concerning the fast dynamics the third level is involved in. So that 
the effective Hamiltonian of the system including the ac-Stark shift, in the dipole and rotating 
wave approximation, can be written as [2D [22] {fi = 1) 

H eff . = tfafaS® + p 2 Sgjj) + tfafaS® + faS$) 

+Ci {S2a 2 + SgP) + C 2 (e^a 2 + e-^^at 2 ) . (6) 

We denote by (3 X and (5 2 the intensity-dependent Stark shifts (3 1 = C 2 /A, an d P2 = Co /A, that 
are due to the virtual transitions to the intermediate relay level and Q = ^ 2 , (A / 0). This 
means that the two three-level trapped ions (one-photon transitions) can be described by by an 
effective two two-level system (in this case two-photon process). A scheme utilizing position- 
dependent ac Stark shifts for doing quantum logic with trapped ions has been presented [23] • It 
has been shown that specific ac Stark shifts can be assigned to the individual ions using a proper 
choice of direction, position, and size, as well as power and frequency of a far-off-resonant laser 
beam. 

The time evolution of the system density operator p(t) can be written as 



where 7 is the phase decoherence rate. Equation (J7j) reduces to the ordinary von Neumann 
equation for the density operator in the limit 7 — > 0. The equation with the similar form has 
been proposed to describe the intrinsic decoherence [27]. Under Markov approximations the 
solution of the master equation can be expressed in terms of Kraus operators [1] as follows 

pit) = jr exp (-iHt) exp (-^H 2 ] p(0) exp (-^H 2 ) exp (iHt) H m 

m=0 m ' V / V / 

= y ^-M m (t)p(0)M^ m (t), (8) 
^— ' ml 

where p(0) is the density operator of the initial state of the system and M m are the Kraus 
operators which completely describe the reduced dynamics of the qubits system, 

M m = H m eM-iHt) exp (~y^ 2 ) • (9) 

Equation (jSJ) can also be written as 

p{t) = exp (-iHtj exp (~y^ 2 ) {e §Mt p(0)} exp ^-^-H 2 j exp (iHt} 



where we have defined the superoperator SmP(0) = Hp(Q)H. We will choose the following 
mixed state of the ions 

p a (0) = cos 2 9 \a, b) (a, b\ + sin 2 9 \b, a) (b, a\ e & A . (10) 

while the initial state of the vibrational mode is in a vacuum state p / (0) = |0> <0| e & F - Then 
the initial state of the system can be written as p(0) = p a (0) (g> p* (0). 

3 Entanglement 

In this paper we take the measure of negative eigenvalues for the partial transposition of the 
density operator. It was proved that the negativity is an entanglement monotone [29j, hence, the 
negativity is a good entanglement measure. According to the Peres and Horodecki's condition 
for separability [30] , a two-qubit state for the given set of parameter values is entangled if and 
only if its partial transpose is negative. The measure of entanglement can be defined in terms 
of the negative eigenvalues of the partial transposition in the following form [31] 

/ P (t) = 2max(0,-A« 9 ) (11) 

where \ff eg is the sum of the negative eigenvalues of the partial transposition of the time- 
dependent reduced atomic density matrix p a , which can be obtained by tracing out the vibra- 
tional mode variables 

p a = Tr f (p(t)), (12) 

In the two qubit system (C 2 <g> C 2 ) it can be shown that the partial transpose of the density 
matrix can have at most one negative eigenvalue [30J. The entanglement measure then ensures 
the scale between and 1 and monotonously increases as entanglement grows. An important 
situation is that, when I p (t) = the two qubits are separable and I p (t) = 1 indicates maximum 
entanglement between the two qubits. In our calculations, we have used the two qubit basis 
\aa), \ab), \ ba) and \bb) to obtain the evolution of the density matrix of the system. 

An interesting question is whether or not the entanglement is affected by the different 
parameters of the present system with the initial state in which one of the qubits is prepared 
in its excited state and the other in the ground state. In particular, we focus on the effect of 
the mixed state parameter 9, the Stark shift parameter (3, (= /3 X ) and the decoherence. As 
expected from the results presented in [32], the analytical solution of equation (jSJ), when we 
set 7 = 0, does not depend on the Stark shift parameter /3 2 , this result can be understood as 
coming from the setting of the initial state of the vibrational mode, which was assumed to be 
in the vacuum state. 
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Figure 1: The evolution of the quantum entanglement I p (t) as a function of the scaled time 
Xt, (A = Ci = C2) an d Stark shift parameter (3. The parameters are 9 = and 7 = 0. 

A numeric evaluation of the entanglement measure leads to the plot in figure 1. We consider 
the initial state of the two ions 9 = 0. In this case, we see that the entanglement is equal to 
zero in a periodic way for a small values of the Stark shift parameter, this period is increased 
with decreasing the parameter f3. It is remarkable to see that with the value of the Stark shift 
parameter, (3 = 20 the entanglement is only zero for the initial period of the interaction time, 
while long lived entanglement is observed as the time goes on. In this case we can say that, 
when the system is allowed to evolve without applying a phase shift (0 = 0), the entanglement 
is a periodic function of time for small values of the Stark shift while long-survival entanglement 
can be obtained for larger values of Stark shift (see figure 1). This is particularly because of 
the nonlinear nature of the coupling in this case (two-photon process) |28j . 




We see from the figure 2a that large values of the Stark shift parameter leads to zero en- 
tanglement. The situation here is quite different from that observed in figure 1, where the 
entanglement exist only for small values of Stark shift while vanishes for all periods of the 
interaction time when (3 > 12. In figure 2b, we pause to touch on certain entanglement fea- 
tures when a mixed state of the qubits is considered as p a (0) = 0.5(|a, 6) (a, b\ + \b, a) (b,a\) 
i.e 9 = 7r/4. It is interesting to see here that the first maximum value of the entanglement is 
observed at earlier time than the previous case. For large values of the Stark shift parameter 
we see that the entanglement has zero value only for a short period of the interaction time 
and then starts to increase. This zero entanglement period is increased when the Stark shift in 
increased further. These properties show that the role played by the Stark shift on the entan- 
glement is essential. Interestingly, when (3 is taken to be non zero, the values of the maximum 
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Figure 2: The evolution of the quantum entanglement I p (t) as a function of the scaled time 
At and different values of Stark shift parameter (3, where, (3 = 2 (solid curve), (3 = 5 (dotted 
curve) and (3 — 15 (dashed curve). The parameters are 7 = and (a) 9 = ~ and (b) 8 = |. 
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entanglement are decreased, indicating that the mixed state setting leads to a decreasing of the 
qubit-qubit entanglement. Generally speaking, because of the influence of mixed state param- 
eter on entanglement, the amplitude of local maxima and minima decrease with increasing the 
deviation of (3 from the unity. However, as f3 takes values close to the unity we return to the 
same behavior in the initial pure state setting i.e p = \a, b) ® (a, b\ . However a slight change in 
/3 therefore, dramatically alters the entanglement. This is remarkable as the entanglement is 
strongly dependent on the initial state, which can be entangled or unentangled. 

We devote the discussion in figure 3 to consider the decoherence parameters effect on the 
entanglement in the presence of Stark shift. We would like to remark that decoherence due to 
normal decay is often said to be the most efficient effect in physics. Which means that, the 
entanglement increases rapidly, then approaches to a minimum value in a periodic manner. In 
this case, the entanglement introduced by the coherent interaction oscillates without dissipation, 
as showed in Fig. 3. Once, the environment has been switched on, i.e., 7 7^ 0, it is very clear 
that the decoherence plays a usual role in destroying the entanglement. Also, from numerical 
results we note that with the increase of the parameter 7, a rapid decrease of the entanglement 
(entanglement sudden death) is shown (in agreement with [1]). 
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Figure 3: The evolution of the quantum entanglement as a function of the scaled time Xt and 
different values of the decoherence parameter 7, where, 7 = 0.01 (solid curve), 7 = 0.1 (dotted 
curve) and 7 = 0.7 (dashed curve). The other parameters are 9 = n/2 and /3 1 //5 2 = 1. 

The remaining task is to identify and compare the results presented above for the entangle- 
ment degree with another accepted entanglement measure such as the concurrence [29]. One, 
possibly not very surprising, principal observation is that the numerical calculations corre- 
sponding to the same parameters, which have been considered in figures 1-3, give nearly the 
same behavior. This means that both the entanglement due to the negativity and concurrence 
measures are qualitatively the same. 
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4 Conclusion 

We have investigated the entanglement in the context of an ensembles of two identical qubits 
(or ions) and negativity as computable measure of the mixed-state entanglement has been used. 
We have treated the more general case where the initial state of the two qubits can be mixed 
taking into account the presence of Stark shift. Through analysis, we find that the extent to 
which that the entanglement vanishes due to Stark shift relies not only on the Stark shift value, 
but also on the initial state setting. When the two ions start from a mixed state, the larger 
the Stark shift is, the faster the entanglement vanishes. For pure quantum states, the complete 
disentanglement occurs for a very short time in a periodic way only for small values of the Stark 
shift. We found that, the entanglement decay due to Stark shift for an initial mixed state is 
similar to the entanglement decay due to the decoherence. Finally, we expect our work will be 
helpful for preserving entanglement in practical experiments. 
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